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The possibility that we are living in a fast expanding and underdensed local bubble has made 
many to debate if dark energy is needed to explain the apparent over-dimming of distant super- 
novae la (SNe la). Recently the Planck team has reported a lower value of Hubble constant and 
a larger matter density than known previously. Here we show that the lower Hubble constant 
is inconsistent with SNe la data, and the local bubble density is too low, unless it is also a global 
property of the universe at the present cosmic time. We suggest a new scenario that the universe 
expands at initially a low, then slightly higher, and finally much higher rate at present, corre- 
sponding to increasing Hubble constant with cosmic time. Therefore these data provide evidence 
for re-accelerating expansion of the universe, deviating from accelerating expansion described 
by the concordant cosmological model, but still not requiring preferred observers. 



Explorations of Type la supernovae (SNe la) have 
proven to be an extremely fruitful enterprize. For ex- 
ample, comparisons between the apparent magnitudes 
of the low and high redshift SNe la have allowed the 
discovery of the accelerating expansion of the universe 
(7, 2). The Hubble constant Hq is the single most im- 
portant cosmological parameter, because it measures the 
expansion rate of the universe, provides the basic infor- 
mation on the age of the universe, and is a key parame- 
ter in determining other cosmological parameters. The 
best measurement of 77o can be made using nearby Type 
la supernovae (SNe la), which are currently the best 
standard candles in cosmology. Recently a 3.3% error 
of Hq is reported by calibrating these standard candles 
with many Cepheid variables in their host galaxies (5), 
which are the best distance indicators of the local uni- 
verse. With such a small error and being independent 
of or at least insensitive to cosmological models, Hq 
can be used directly as a pre-determined parameter in 
determining other cosmological parameters (4). Being 
excellent standard candles, SNe la can also be used to 



calibrate other candles to extend the distance ladder to 
much higher redshift, such as gamma-ray bursts (5) and 
super-Eddington accretion black holes (6). 

However, the best measurement of Hq was made us- 
ing eight nearby SNe la with redshift z from 0.0043 
to 0.0072, yielding only a local Hubble constant Hq q 
(at z ~ 0), which may not necessarily be the same as 
the Hubble constant at a much higher redshift, Ho, z (at 
z 0; note that it is not the same as the commonly 
used Hubble parameter H{z)). Indeed, just before the 
discovery of the accelerating expansion of the universe 
(7, 2), evidence was found that TTo.o > Ho,z by about 
6%, where the boundary is around Dl ~ 70 /i _1 Mpc 
and h = H o /100 km s" 1 Mpc" 1 (7). This suggests 
that we are living within a small local Hubble bubble of 
radius about 70 /i _1 Mpc, which expands slightly faster 
than the outside universe. This means that we are mov- 
ing away with respect to distant SNe la faster than the 
global Hubble expansion and thus the distant SNe la 
should look dimmer than viewing only within the Hub- 
ble flow. This thus has led many to debate if the ac- 
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celerating expansion of the universe is simply an mi- 
rage of this local Hubble bubble (8-16), since the over- 
dimming of distant SNe la is what led to the initial dis- 
covery of the accelerating expansion of the universe, 
which has been explained as due to the yet unknown 
dark energy (1, 2). 

Similar to the previous work (7), in Figure 1 we 
show Ho measured from the currently best available 
SNe la data, i.e. the eight SNe la used to obtain 77o,o 
and the Union 2.1 compilation (77) used here to mea- 
sure Ho z . The data for the eight SNe la are listed in 
Table 1; ho (with standard error <jh ) is calculated using 
Equation (4) and data in Table 3 of Ref. (3). ho for other 
SNe la is calculated using the standard method, e.g., 
Equation (12) or (13) in Ref. (18) with the same choices 
of parameters. We limit the SNe la with z < 0.04 (with 
a median redshift of z\ = 0.025), in order to avoid 
any coupling with cosmological parameters; in fact it 
is already safe to choose z < 0.1 (18). We obtain 
Vo = 0.738 ± 0.0155 and h , Zl = 0.704 ± 0.0051, 
and i^o.o > Ho, Zl at 96.4% confidence level (see the 
figure caption for details). Following the Hubble's law 
V r = HqR, in which V T is the receding velocity of a 
distant galaxy and R is its radial distance from the ob- 
server, the above result indicates that our local universe 
is expanding at a greater rate than the distant universe. 
This confirm the previous result (7) with the most up- 
dated and best available data. 

A striking new result just released by the Planck 
team is that M = 0.307 ± 0.019 and h , Z2 = 0.679 ± 
0.015 and (z 2 ~ 1100) (79), the latter of which is even 
lower than ho, Zl , as also shown in Figure 1. This result 
is actually consistent with some earlier reports (20, 21). 
In Figure 2, we examine the consistency of the three dif- 
ferent Hubble constant, i.e., ho,o, ho. Zl and ho, Z2 against 
the Union 2.1 SNe la data at both low and high redshift, 
within the framework of the ACDM model with differ- 
ent values of Hm- As expected, only ho = 0.704 is 
consistent with data at low redshift (z\ = 0.025), inde- 
pendent of S1m- At high redshift (z 2 = 0.740), again 
only ho = 0.704 is consistent with data unless S7m de- 
viates significantly from the well-measured and widely 
accepted value of around 0.3; actually the high redshift 
SNe la data favors a slightly lower fijyr* consistent with 
the definition of SIm = 8nGpo /H 2 f° r a constant rhoo- 
It is thus very unlikely that SNe la data can be recon- 
ciled with either the higher or lower ho measured in the 
local bubble or with the cosmic microwave background 
data with Planck, respectively. 



From each listed Ho.o (measured in the local bub- 
ble) and its error in Table 1, the pure statistical error 
of ho,o should be 0.009, much smaller than the error of 
0.0155 determined from the variance of the eight data 
points. This means that the probability that the data do 
not contain additional fluctuations is less than 0.68%. It 
has been known that the peculiar motions of the SNe la 
hosts may cause such fluctuations beyond the measure- 
ment statistical errors (22, 23). For a peculiar velocity 
(Vlos) along the line of sight (LOS), the deviation to its 
cosmological redshift (z) is Az = V(l + z)/c, where c 
is the speed of light. For Vlos = 100 km s _1 at z -C 1, 
Az = 0.0003, comparable to redshift z = 0.0043 to 

0. 0072 of the eight SNe la. Since H = cz/D L when 
z = 0, a non-negligible deviation to 77o may therefore 
be produced, if one takes the common practice of using 
its luminosity distance D^to measure 77o- 

The additional fluctuations caused by the random 
peculiar motions of the SNe la hosts can be found 
from a 2 ho p = {hi) - a\ Q = 0.0352 2 , where (h 2 ) = 
J2(ho,i~h ) 2 /(n-l)anda ho = E^h 0li /n = 0.0262 
(n = 8); in fact a ho ~ a hoi . Clearly a hop > a ho , 

1. e., the average fluctuation to ho caused by the putative 
random peculiar motions of SNe la hosts is larger than 
the measurement errors in ho- This means that for the 
first time it is possible to measure the peculiar motion of 
each individual SNe la host with the extremely accurate 
measurement of ho : i of each SNe la. 

We then compare the measured |77o,i — Hq\ with the 
expected deviations caused by different LOS peculiar 
velocities at low redshift in Figure 2; here 77o = 77o,o- 
The data are consistent with Vlos ~ 100 km s _1 . In 
Figure 3, we show the positions of the eight SNe la in 
equatorial coordinates and their LOS peculiar velocities 
from the Hubble flow (all these data are listed in Ta- 
ble 1); the detected peculiar motions do not show any 
significantly coordinated pattern (albeit with small num- 
ber statistics) and thus are consistent with random mo- 
tions with respect to the Hubble flow. For random pecu- 
liar motions, (Vp oc ) (the variance of its space velocity 
Vp ec ) is related to (Ulos) (the variance that of the its 
LOS velocity Vlos) by (Vp 2 cc ) = ^(V 2 oc ), where 
n is the number of SNe la, and x follows a gamma- 
distribution, 

p( X ;a,l3) = -^-x a - 1 e^ x , (1) 
T(a) 

where T(a) is the commonly known T-function, a and 
j3 are two parameters: a = ^(n — 1) and (3 = fto + 
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Table 1 : Redshift, coordinates (J2000.0), Hubble constant and LOS peculiar velocities of the eight SNe la. 



Name 


redshift 


Right ascension 


Declination 


Hubble Constant 


LOS pec. vel. 






(hh:mm:ss) 


(dd:mm:ss) 


H (kms- 1 Mpc" 1 ) 


Vlos (kms" 1 ) 


1981b 


0.0072 


12:34:29.57 


+02:11:59.3 


70.37 (2.26) 


-100.8 (66.7) 


1990n 


0.0043 


12:42:56.68 


+13:15:23.4 


74.81 (2.75) 


17.8 (48.3) 


1994ae 


0.0045 


11:56:25.87 


+55:07:43.2 


76.38 (2.46) 


47.3 (45.1) 


1998aq 


0.0055 


10:47:01.94 


+17:16:30.8 


70.89 (2.61) 


-65.2 (58.6) 


1995al 


0.0059 


09:50:55.97 


+33:33:09.4 


76.87 (2.83) 


74.1 (68.1) 


2002fk 


0.0070 


03:22:05.71 


-15:24:03.2 


68.48 (2.83) 


-151.9 (81.1) 


2007af 


0.0070 


12:01:52.80 


-18:58:21.7 


81.43 (2.62) 


218.2(75.0) 


2007sr 


0.0063 


14:22:21.03 


-00:23:37.6 


70.50 (2.59) 


-84.8 (66.8) 



l/2\/2(n - 1), where /3 = 0.68269. For n = 8 here, 
the 90% lower limit of a; is £ > 0.462. From the eight 
SNe la, we have y/Whos) = 91.2 km s -1 (the mea- 
surement errors of Vlos have been subtracted), which 
gives the expectation, most-likely and 90% upper limit 
values of y/(V£ os ) as 163.4, 266.6 and 355.0 km s"\ 
respectively. 

Peculiar motions of galaxies can also be measured 
indirectly with the two-point correlation function to 
quantify one of the redshift space distortions, i.e., the 
apparent elongation of the clustering of galaxies in red- 
shift space along LOS, the so-called "fingers of God" 
effect caused by peculiar motions of galaxies. This was 
first tentatively observed with Vp ec ~ 200 — 600 km s _1 
for nearby galaxies (24) and later convincingly demon- 
strated with Vp cc ps 400 km s _1 in the 2dF redshift sur- 
vey at a median redshift of z = 0.11 (25). The 2dF 
result is consistent with £Im ~ 0.3. Therefore the ob- 
served peculiar velocities of the eight SNe la suggest 
that the expectation, most-likely and 90% upper limit 
values of the local matter density f2M,Loc are 0.04, 0.11 
and 0.20, respectively, according to the Cosmic Virial 
Theorem (24), i.e., y/(V£ os ) ps 800S1^ 2 km s" 1 . Our 
local expanding bubble is thus most-likely underdensed 
significantly, compared to the average matter density of 
the universe. 

If the observed ho,o = 0.738 is caused by the under- 
densed local bubble embedded in a Hubble expansion 
with ho = 0.704, then 51m. Loc should be only about 
10% lower than the global Qm, and thus significantly 
larger than that measured above with the observed pe- 
culiar velocities of the eight SNe la. In other words, if 
the most recent universe expands at the same rate as the 



early universe, our local Hubble constant should be even 
much larger than ho,o = 0.738, due to the significant 
under-density in the local bubble. A way out is that the 
local bubble is actually the global property of the uni- 
verse; observationally it becomes "local" because only 
a small volume of the most recent universe can be ob- 
served by any observer. In other words, an observer lo- 
cated anywhere in the universe at z ~ 0, with respect to 
the cosmic microwave background at z ~ 1100, should 
also observe the same local bubble. This is good news, 
since it naturally avoids the philosophical crisis if we are 
living in a specially chosen place in the universe, i.e., the 
center of the universe where the matter density is much 
lower. 

Therefore the SNe la and Planck data support a new 
scenario that the universe expands initially at a low rate 
(at z ~ 1100), then at a slightly higher rate (at z < 1), 
and finally at a much higher rate at present (at z ~ 0). 
We call this re-accelerating expansion of the universe, 
to distinguish it from the well-known accelerating ex- 
pansion of the universe described by the ACDM model, 
dominated by dark energy and with a constant Hubble 
constant (1, 2). Mathematically this model of universe 
can be described with the well-known the Lemaitre- 
Tolman-Bondi (LTB) model with A, in which the Hub- 
ble constant, and perhaps also other cosmological pa- 
rameters, are functions of cosmic time only; however, 
it is not clear what drives the time varying Hubble con- 
stant and thus the re-accelerating expansion of the uni- 
verse. Nevertheless, it has been shown that the Union 
2.1 SNe la data also agree with the LTB model, but 
cannot distinguish between the LTB model and ACDM 
models alone (26). Finally we stress that this new sce- 
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nario can still maintain the Copernicus principle in the 
way that there is no preferred observer in the universe. 
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Figure 1: The Hubble constant ho measured with each SNe la within a luminosity distance Dl of 250 h~ 1 Mpc. 
Left panel: The red crosses are the eight best SNe la (3) used to measure the local Hq with Dl < 50 h~ 1 Mpc, 
giving an average local ho,o = 0.738 ± 0.0155 marked as the thick solid red line. The blue crosses are the Union 
2.1 SNe la (77) at D L > 80 ft^Mpc, yielding an average cosmological h , Zl = 0.704 ± 0.0051 as the thick 
solid blue line {z\ = 0.025 is the median redshift of these SNe la marked by the blue crosses). Seven of the 
eight SNe la (red crosses) have ho > ho tZl , indicating that the probability that the eight SNe la are drawn from 
the same population of the other SNe la (blue crosses) is less than 3.6%. Right panel: Histogram of the blue 
crosses in the left panel. The filled red and blue areas are the ler error regions of /io,o an d ^o,zi respectively; their 
errors are calculated from the variance of each sample, and are significantly larger than that calculated from error 
propagation using the measurement errors of all data points (see text for details). The large error in ho,o is due to 
its very small sample size of only eight data points and additional fluctuations caused by the peculiar motions of 
their hosts (see text for details). hop an d ho, Zl are different at 2.1c level with respect to their joint error bar, i.e., 
the probability that they are consistent with each other is less than 3.6%. For comparison, the just released Planck 
result h Q . Z2 = 0.679 ± 0.015 is also marked by the filled magenta area (z 2 ~ 1100). 
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Figure 2: Residuals of the distance modules (/i) of the Union 2.1 SNe la against the prediction of the ACDM 
cosmological model with different combinations of Qm and and Hq, J7m + = 1 is always assumed. The black 
thick vertical dashed lines indicate A/i = 0. Two groups of SNe la are chosen here: low redshift of z < 0.04 with 
a median redshift of 0.025, and high redshift of z > 0.5 with a median redshift of 0.740. The combination of the 
absolute value of A/i (the center of the distribution of A/i) and the magnitude of cr^ (the standard deviation of 
each Gaussian fit, also labelled in each panel), indicates how well the model describes the data. Each filled area 
marks the 3<r error range of A/i with cr^ = oa^/Vh — 1, where n is the total number of data points. The four 
panels plotted in magenta color are consistent with data: ho = 0.704 for z = 0.025, independent of Qm; also 
ho = 0.704 for zi = 0.740, unless Hm deviates significantly from 0.3. 
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Figure 3: Deviations of the measured Hq with individual SNe la, caused by the peculiar motions of SNe la hosts 
at different redshift. The lines are the expected deviations for different peculiar velocities along the LOS. The data 
points with error bars are the deviations of the eight SNe la from /io,o = 0.738. The data are consistent with the 
LOS peculiar velocity of about 100 km s _1 . 
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Figure 4: Positions of the eight SNe la in equatorial coordinates and their LOS peculiar velocities from the Hubble 
flow. Negative (marked as blue diamonds) and positive (marked as red triangles) velocities mean that their hosts 
are moving towards and away from the observer within the Hubble flow with /lo.o = 0.738, respectively. The sizes 
of the signs are proportional to their LOS velocity deviations. 
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